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We investigate models for the photoionization of the widespread diffuse ionized gas in galaxies. 

In particular we address the long standing question of the penetration of Lyman continuum pho- 
^Nj tons from sources close to the galactic midplane to large heights in the galactic halo. We find that 

recent hydrodynamical simulations of a supernova-driven interstellar medium have low density 
.^ paths and voids that allow for ionizing photons from midplane OB stars to reach and ionize gas 
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many kiloparsecs above the midplane. We find ionizing fluxes throughout our simulation grids 
are larger than predicted by one dimensional slab models, thus allowing for photoionization by O 
stars of low altitude neutral clouds in the Galaxy that are also detected in Ha. In previous studies 
of such clouds the photoionization scenario had been rejected and the Ha had been attributed 
to enhanced cosmic ray ionization or scattered light from midplane H li regions. We do find 
that the emission measure distributions in our simulations are wider than those derived from Ha 
Ci observations in the Milky Way. In addition, the horizontally averaged height dependence of the 

gas density in the hydrodynamical models is lower than inferred in the Galaxy. These discrep- 
^~^ ancies are likely due to the absence of magnetic fields in the hydrodynamic simulations and we 

JZi discuss how magnetohydrodynamic effects may reconcile models and observations. Nevertheless, 

■^ we anticipate that the inclusion of magnetic fields in the dynamical simulations will not alter our 

CO primary finding that midplane OB stars are capable of producing high altitude diffuse ionized 

gas in a realistic three-dimensional interstellar medium. 
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1. Introduction 

Diffuse ionized gas (DIG, also commonly referred to as the Warm Ionized Medium in the Milky Way) 



is observed along all lines of sight in the Milky Way (Haffner et al. 2003). Observations of Ha emission 
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(Haffner et al. 1999), pulsar dispersion measures (Reynolds 1989 Gomez et al. 2001 Gaensler et al. 2008 



Savage fc Wakker|2009| , and [Al III] absorption (Savage fc Wa kker|2009 1 indicate the DIG in the Galaxy has 
a scale height of 1 — 1.5 kpc and an H+ column density about 1/3 that of the neutral hydrogen (Reynolds 
1991). Widespread DIG is a common feature in spiral galaxies. It is detected to large heights in edge-on 



galaxies ( Rand|[l998 ) and in the inter-arm regions of face-on galaxies well away from traditional H II regions 
(Ferguson et al. 1996). For a summary of the physical properties (ionization state, electron temperature, 



scaleheight) of the DIG in the Milky Way and other galaxies see the recent review by Haffner et al. (2009). 



The power requirements of the DIG indicate that the only reasonable source of ionization is from O 
stars (Reynolds 1990a|b ). However, in a smoothly distributed interstellar medium (ISM) with a hydrogen 
number density of 1 cm~'^, an O star with an ionizing luminosity Q = 5 x 10*^ s~^ can only ionize a volume 
of radius 53 pc (assuming a hydrogen recombination coefficient as = 2.59 x 10"^"^ cm^ s~^). Photoionization 



models including O stars in a vertically stratified ISM (e. g. [Miller fc C"ox||1993[ [Dove fc Shun||1994[ [Wood 
fc Mathis|2004 Ciardi et al.|2002 ) show that ionizing photons can reach to large distances perpendicular to 
the midplane, but in general require a lower Hi density than inferred from observations. Clumping in a three 
dimensional (3D) ISM naturally provides low density paths for ionizing photons to reach large heights above 
the midplane. Indeed, models that adopt an average ISM density re-arranged into a 3D fractal structure 
readily allow photons to penetrate from the midplane to large heights (e. g. Ciardi et al.|2002" Haffner et al. 
2009). Turbulence is perhaps the most plausible mechanism for producing such clumping and the DIG is 



well known to be turbulent ( Armstrong et al!]|1995 Benjamin||1999 Hill et aL||2008 Chepurnov & Lazarian 
2010). The concept of 3D structures allowing deeper penetration of radiation compared to smooth media 



is not new or unique to photoionization studies. Previous theoretical work has demonstrated the increased 



mean free paths for non-ionizing photons in 3D dust structures in the ISM and molecular clouds (Boisse 
1990), stellar winds (,Shaviv^^2000j) , and exoplanetary atmospheres (Hood et al.||2008). 



In this paper we extend our investigation of O star photoionization modeling of a 3D ISM to study 



photoionization of 3D hydrodynamic simulations of supernova-driven turbulent ISM (Joung & Mac Low 



2006 Joung et al. 20091. In § 2 we briefly describe the hydrodynamical and photoionization simulations; 



§ 3 presents the results of our photoionization simulations showing the resulting distributions of ionized and 
neutral gas; in § 4 we discuss discrepancies between our models and observations of the Galactic DIG and 
the likely role of magnetic fields not present in the current simulations. 



Models 



2.1. Hydrodynamic Simulations 

In this paper we investigate the ionization structure of a 3D density field from a hydrodynamic simulation 
of a supernova-driven turbulent ISM ( Joung fc Mac Low||2006 Joung et al.||2009 |. The simulations use the 
FLASH v2.4 code (Fryxell et al. 20001. We summarize the essential features of the models here. The adaptive 



mesh simulation has a maximum resolution of just under 2 pc in a region with \z\ < 250 pc, and progressively 
coarser resolution at larger altitudes. The full-resolution simulation box has dimensions 0.5 kpc in x and y 
and extends to ±5 kpc in z above and below the midplane. A fixed gravitational potential from |Kuijken] 



& Gilmore ( 1989 1 was used. The model employed here, with the Galactic supernova rate, has Type II and 
Type I supernova rates of 27.4 and 6.58 Myr^^ kpc^^ and scale heights of 90 and 325 pc, respectively. Of 
the Type II supernovae, 3/5 are concentrated spatially and temporally to simulate superbubbles and their 
heights are confined to \z\ < 90 pc. The total gas surface mass density is 1.87 x 10^ M0 kpc~ . 



A diffuse heating term representing photoelectric heating ( Wolfire et al. 1995 ) of F 



8.50x10-26 ergs-i 

and radiative coohng appropriate for an optically thin, solar metallicity plasma are included. The hydrody- 
namical simulation does not track ionization of the gas, so photoionization heating is not explicitly included. 
In the DIG, photoionization by a dilute radiation field is likely the dominant heating mechanism, although 
there may be an additional mechanism that is important at low densities, n < 0.1 cm~^ (Reynolds et al 



1999 Wood & Mathis 2004). However, this has a small impact on the dynamics of the ISM because the 



thermal pressure is small compared to the turbulent pressure. In addition, Selkowitz & Blackman (2007) 



have shown that magnetic fields may provide additional heating in the DIG. This is very relevant especially 
given our conclusions on the need for the inclusion of magnetic fields in the dynamical simulations to explain 
the observed emission measure distributions in the Galactic DIG. 

In their hydrodynamic simulations of a vertically stratified ISM driven by individual supernova explo- 
sions, Joung & Mac Low (2006) find that the density power spectrum peaks at ~ 20 pc and that the most 
energy-containing scales lie at 20-40 pc (see their Fig. 8). On the other hand, they find no single effective 
driving scale; energy injection occurs over a broad range of scales, with greater than 90% of the total ki- 
netic energy contained in wavelengths shortward of 200 pc. Spatially correlated supernovae explosions that 
break out of the gaseous disk (i.e., superbubble blow outs) and high-altitude Type la supernovae are the 
main drivers of turbulent motions of the halo gas. These structures have larger characteristic scales (several 
hundred parsecs) than individual supernova remnants in the disk, due to the higher amount of energy and 



lower background gas density, respectively; see Fig. 2 of Joung & Mac Low (2006). 



For our photoionization simulations we consider a subset of this grid extending to ±2 kpc in z. Memory 
requirements of our photoionization code on current desktop computers required us to re-bin the density grid 
by a factor of four giving a resolution of 7.8 pc per grid cell. The computationally straightforward approach 
of running the radiative transfer code in a static realization of the hydrodynamical code is justified by the 
recombination time scale of ~ 3 Myr, much shorter than the dynamical timescale of ~ 100 Myr in an 8000 K 
hydrogen gas with a total number density of 1 cm^^ and rig = 0.03 cm"*^ ( Spitzer|l978 ). 

Density, velocity, and temperature slices of the simulation are shown in Joung et al. (2009). In Fig. 1 



we show the horizontally averaged density as a function of height compared to a model for the hydrogen 



density in the ISM. The model comprises a Dickey-Lockman (Dickey & Lockman 1990) distribution for 



neutral hydrogen with a low density, vertically extended distribution for ionized hydrogen as follows: 
1 
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where the number density is in units of cm^'^ and the heights, z are in pc. The first three terms represent 
the density of neutral hydrogen and the fourth term represents ionized hydrogen. The mean vertical column 
density (within the range \z\ < 2 kpc) for the simulations is in good agreement with that derived from 
equation 1 with a value of ~ 7 x 10^° cm^^. However, the density structure of the dynamical simulation is 
more centrally peaked than the Galactic model described by equation 1. As we discuss later, this is most 
likely due to the absence of magnetic fields in the current dynamical simulations. 



2.2. Photoionization Simulations 



The primary focus of this paper is the ionization structure of hydrogen and the penetration of Lyman 
continuum photons to large heights above the Galactic midplane. Therefore for the photoionization simu- 



lations we use the 3D hydrogen-only Monte Carlo photoionization code described in Wood & Loeb (20001 



Because neutral hydrogen is opaque to Lyman continuum radiation whereas the optical depth of nearly 
fully ionized hydrogen is much lower, we use an iterative process to establish a self-consistent distribution 
of neutral and ionized hydrogen. This code follows the propagation of ionizing photons within a 3D linear 
cartesian density grid and accurately computes the 3D ionization structure of hydrogen for an assumed 
uniform temperature. The code includes photoionization arising from both direct stellar photons and also 
the diffuse radiation field produced by recombinations to the ground level of hydrogen. We have extended 



the original Wood & Loeb ( 2000 ) code to include the effects of dust scattering and absorption of ionizing 
photons. We adopt an opacity of 1.3 x 10~^^cm~^H~^ appropriate for silicate dust at Lyman continuum 



wavelengths as discussed in Mathis & Wood ( 2005 1 , and assume albedo and scattering phase function asym- 
metry parameters a = 0.5, g = 0.5. While dust opacity may influence the ionization structure of high density 
H II regions, its effects in the low density DIG that we are investigating are not important. We find that 
typically less than 1% of ionizing photons are absorbed by dust in our simulations. The simulated emission 
measure distributions in § 3 below do not include regions of the grid with \z\ < 150 pc, thereby further 
mitigating the effects of dust in our analysis. For our current study we are not investigating emission line 
ratios, so the single temperature hydrogen-only code is sufficient. It has the advantage of being faster and 
requiring less memory than our code that computes the ionization structure of many more elements and also 
the electron temperature ( Wood et al.||2004 ). 

The sources of ionizing radiation are point sources within the simulation grid representing single stars or 



clusters of OB stars. We follow the approach of Ciardi et al. (2002) and randomly choose x and y locations of 



the sources uniformly within the grid. The z coordinate of each source is randomly chosen from a Gaussian 
distribution with a scale height of 63 pc ( Maiz-Apellaniz|[2001 ) . 



The catalogue of 424 O stars within 2.5 kpc of the Sun (Garmany et al. 1982) indicates a stellar 
surface density of 24 stars kpc~^. The Lyman continuum luminosity from these stars was estimated to be 
Q — 7 X 10^^ s~^ (Vacca et al. [19961. The density field from the hydrodynamic simulation represents a slice 



of the ISM with a horizontal width of 500 pc, so adopting the estimated stellar surface density, we would 
expect around six sources within our grid providing an ionizing luminosity of around 10^*^ s~^. A fraction of 
these ionizing photons will produce the DIG, with the rest producing local H ii regions or escaping from the 
simulation entirely. We do not consider the formation of high density H ii regions around sources, so our 
models are really considering the ionizing photons that have escaped from their local environment into the 
lower density diffuse ISM. In our simulations we therefore select random locations for six ionizing sources. 
We uniformly distribute the ionizing luminosity among the six sources and investigate total luminosities 
(i.e., ionizing photons escaping from traditional H ii regions) in the range 0.1 < Q^g < 10, where Q49 = 
Q/(10'*^ s~^), covering a range of scenarios allowing for a lower and higher Lyman continuum budget than 



estimated by Vacca et al. ( 1996 ) 



Repeating boundary conditions are adopted for the Monte Carlo radiation transfer whereby ionizing 
photons that exit an x or y face of the grid are re-injected on the opposite face of the grid with their 
direction of travel unaltered. This simulates the more realistic situation within the Galaxy of ionizing 
photons that originate from distant sources. Therefore ionizing photons can only exit through the top or 
bottom z faces of the grid or if they are degraded to a non- ionizing photon or absorbed by dust. 



3. Results 

We investigated a large number of models for different source locations and total ionizing luminosities. 
The position of the sources within the density grid and the total ionizing luminosity affect the vertical extent 
to which gas in the hydro simulation can be ionized. For some realizations of the randomly sampled source 
locations, sources may be placed many tens of pc above the midplane, and thus above the densest regions in 
the hydrodynamic simulation. Such a positioning of sources results in more of the gas at large heights being 
ionized compared to simulations where the sources are confined to midplane regions. Similar results were 



found in the ionization simulations of Ciardi et al. (2002) who investigated the escape of ionizing photons 



from vertically stratified smooth and 3D (fractal) ISM density distributions. 

The single parameter that most strongly influenced our results and determined the vertical extent to 
which the gas could be ionized was the total ionizing luminosity, Qtg. The number and location of ionizing 
sources did not have a big effect - if a source was randomly placed within a very dense region of the grid and 
its ionizing photons were trapped, there were always other sources in less dense regions that could provide 
the photons to ionize gas at large heights. In the following sections we quantitatively present the results of 
our photoionization simulations. 



3.1. Distribution of neutral and ionized hydrogen 

Figure 2 shows the average neutral fraction (np'^J/np^ + H+]) as a function of height for different 
ionizing luminosities. At low values of Q49 < 1, the grid becomes fully neutral for \z\ > 800 pc. We find 
that the grid is almost fully ionized at large heights (|z| > 0.5 kpc) for total ionizing luminosities Q4g > 2. 
Figure 3 displays the emission measure (EM = j n^ds) and H I column density maps for edge-on viewing of 
the grid, showing the development of extended ionized volumes as the ionizing luminosity increases. 

To maintain the ionization of the Galactic DIG in the Solar neighborhood (within 2-3 kpc) requires a 
Lyman continuum flux in the midplane of around 4 x 10^cm~^s~^ (Reynolds 1990al, about 12% of that 



available from OB stars (Abbott 1982 Vacca et al. 1996). The required ionizing luminosity, both for the 



Galactic DIG and for our simulations, is lower than that estimated to be available from OB stars in the 
Galaxy and represents Lyman continuum photons that leak out of traditional (high density) H 11 regions 
to ionize the DIG. As discussed in § 2.2, assuming the Vacca et al. (1996) value for the ionizing flux in the 
solar neighborhood (3.4 x lO^cm^^s"^), the total flux available to ionize our grid is around 10^° s~^. In 
our simulations the gas is ionized at large heights when Q^q > 2, slightly higher than the 12% suggested by 



Reynolds (1990a I 



Figure 4 shows the average neutral and ionized hydrogen densities as a function of height for the 
simulation with Q49 = 2. As described above, the mean densities at large \z\ are lower than the model of 
equation 1 for hydrogen density in the Galaxy. Therefore it is not surprising that the mean density of neutral 
hydrogen for \z\ > 200 pc in the ionization simulations is lower than the Dickey-Lockman distribution. The 
mean electron density from the simulations is also lower than derived for the DIG in the Galaxy, but does 
follow an exponential distribution at large \z\ with a scale height of around 500 pc. This is smaller than the 
scale height of the DIG in our Galaxy. 

For comparison, we investigated ionization of a smooth density grid with the hydrogen density described 
by equation 1. In the simulations we adopted an ionizing luminosity Q^g — 2 and multiplied the first three 
terms in equation 1 by a constant factor until we achieved ionization of the gas at large heights and a total 



escape fraction of ionizing photons of around 10% through the top and bottom of our grid (z — ±2 kpc). 
To achieve this we find that the three components representing the "neutral hydrogen" in equation 1 must 



be multiplied by about 1/3. This is consistent with the ionization study of Miller & Cox ( 1993 ) who found 
that if they used the known locations of O stars in the solar neighborhood then the hydrogen density within 
the smooth component must be lower than the Dickey-Lockman distribution to allow gas at large heights 
to be ionized. Also see the results of 3D simulations by Ciardi et al. ( 2002[ ) and the discussion in Haffner 
et al. (2009). For O stars to provide the ionization of the DIG requires a 3D distribution of gas in which 



the density of any smooth component should be at most one third of the Dickey-Lockman distribution for 
neutral hydrogen. To make up the remainder of the observed neutral hydrogen requires higher density 
"clouds" ( Miller fc Cox||199"3 l or denser regions such as produced in a turbulent medium. 



We further investigated the survival of neutral clouds at large heights above the midplane and performed 
simulations to study the ionization of a 50 pc radius spherical blob positioned at z = 1.5 kpc in the grid. 
Figure 5 shows a slice through the center of the ionization grid containing a blob of various densities showing 
the transition from almost fully neutral to fully ionized with decreasing blob density. In our simulations with 
Q49 = 2, a blob with n > 0.1 cm~'^ remains neutral, except for a very thin ionized skin, not resolvable in our 
current simulations. Lower density blobs exhibit a thicker ionized skin and become neutral deeper into the 
blob. The ionized skin is mostly directed towards the midplane, but does extend all around the blob due 
to photons from the diffuse ionizing radiation field. The number density of neutral hydrogen blobs in our 
simulations agrees well with density estimates (0.05 < n < 0.5) for high velocity clouds (e. g. Wakker & van 
Woerden|p97| [Wakker et aL][2008l [Hill et al.|[2009| . 



3.2. Emission measure distribution 



The results of the previous section demonstrate that the hydrodynamical simulations produce density 
structures that readily allow for ionizing photons to percolate to high altitudes and produce extended layers 
of ionized gas. In Figure 6 we present probability distributions of emission measure, EM = Jn^ds, for 
face-on and edge-on views of the simulation grid. In forming these EM distributions we have removed the 
contributions from grid cells in the simulation with \z\ < 150 pc, which contain the densest regions of 
the simulation and in reality should also include high density H II regions. In forming the EM from the 
simulations we use the computed electron density and not an EM derived from simulated Ha observations. 
An Ha derived EM would require an accurate temperature calculation incorporating heating and cooling 
from a radiation-magnetohydrodynamical simulation, beyond the scope of the current paper. 

For comparison we also show EM distributions in the Galaxy derived from subsets of the Wisconsin 
H-Alpha Mapper (WHAM) Northern Sky Survey ( [Haffner et~aLl|2003[ ). The two EM sin|6| distributions 
displayed are of high-latitude (|6| > 10°) sighthnes that sample 1) all DIG gas north oi S — —30° and 2) the 
Perseus Arm, which is conveniently separated out in the velocity range —100 km s~^ < I^sr ^ —25 km s~^ 



Haffner et al.[[1999). In both cases, we have excluded the large classical H II regions identified by Hill et al. 



( 2008 ) to leave a sample consisting primarily of DIG emission. 



Figure 6 shows that the simulations produce EM distributions that are broader than observed for the 
Galactic DIG. In making these comparisons we compare the face-on view of the simulation with the EM 
distribution for all Galactic sightlines b > 10°, and comparing the edge-on view with the EM distribution 
for the Perseus Arm. The very extended low EM tail seen in the edge-on view of the simulation grid will be 
diminished somewhat in Ha derived EMs since hot (T > 10^ K), low density (n < 10~^cm"^) gas at high 



latitudes (see Joung et al.||2009" ) will have a very low Ha emissivity. 



The standard deviation of log EM in the face-on view of the photoionization simulation is 0.56. This 
is much wider than the observed value of CTEMsiniji = 0.190 for the DIG (Hill et al. 2008). Indeed, in the 



isothermal MHD models presented by Kowal et al. ( 2007 1 and Hill et al. ( 2008 ) , only the highly supersonic 



models produce such a wide distribution of emission measures. Turbulence in the DIG is likely mildly 
supersonic, with a Mach number of 1 — 3 ( Padoan et al.|1997 Hill et al.|2008 Berkhuijsen fc Fletcher|2008 1 . 
In the hydrodynamical models presented here, the warm gas, defined as 5000 K < T < 20000 K, has a mean 
sonic Mach number of Mg = 1.2 — 1.7, as shown in Figure [7j 

The breadth of the EM distributions indicates that the hydrodynamical simulations show much larger 
contrasts between low and high density gas than present in the Galaxy. This result, combined with the 
average hydrogen density at large heights being lower than derived in the Galaxy (Fig. 1), leads us to 
speculate that the effects of magnetic fields (not present in the simulation) could reconcile the model with 
observations. Including magnetic fields in the hydrodynamical simulations will provide two main effects: 
a higher signal speed, and thus a lower effective Mach number, which will result in less extreme density 
variations (as seen in the simulations by de Avillez fc Breitschwerdt||2005 1 and additional pressure support 
that will levitate and maintain a higher density with height. We anticipate that these combined effects 
will result in narrower EM distributions and higher densities with height, bringing our simulations closer to 
observationally derived values for the EM and height dependence of the density of the neutral and ionized 
gas. 



4. Comparison to previous work 

The escape fractions of ionizing photons from our simulations are around 10% for Q = 2 x 10^^ s"^ and 



increase for larger ionizing luminosities, entirely consistent with the study by 



ours and other 3D simulations differ from the ID model of Bland-Hawthorn & Maloney ( 1999 1 that predicted 



Ciardi et al. 



(2002). However, 



very small ionizing fluxes at small distances above the galactic midplane (see Fig. 6 in 'Putma n et al.|[2003 1 . 
This difference is not surprising for several reasons: the Bland-Hawthorn & Maloney ( 1999) model considered 



all ionizing sources to be located in the Galactic midplane and to suffer extinction due to a uniform density 
slab of dust. This is entirely reasonable for their study of ionization of high velocity clouds at very large 
distances from the Galaxy, and for HVCs their model predicts appropriate escape fractions and ionizing 
fluxes. However, for studying photoionization of clouds at small distances from the Galactic plane their 
model predicts a vanishingly small fiux, which is clearly not correct for a 3D simulation. Considering these 
differences between ID and 3D models, we believe it is inappropriate to reject photoionization of the Galactic 

Indeed, the observed 



cirrus as a source for at least some of the observed Ha emission ( Mattila et al 



2007h 



Ha intensities of typically 2 Rayleighs from the clouds studied by Mattila et al. ( 2007 1 require ionizing fluxes 
of around 5 x 10^ photonscm"^ s~^, well within the range suggested for the Solar neighborhood by Vacca 
et al. (19961. For a slab cloud with a uniform density of lOcm"'^ normally illuminated on one side and 
assuming an ionized gas temperature of 8000 K, the ionized skin would be ~ 0.05 pc thick and contribute to 
the observed Ha intensity in addition to any scattered light component. Allowing for a photoionized skin on 
the surface of these clouds means that the observed Ha emission need not be attributed entirely to cosmic 
rays ( del Burgo fc Gambresy|[2006 1 or scattered light from midplane H ii regions (Mattila et al.|[2007). 



Recent models of the DIG in M51 (Seon 2009) suggested that O stars were not capable of producing 



ionization at large distances from their immediate surroundings. This conclusion was based on an ionization 



model of a smoothly distributed ISM (uniform slab or an exponential disk) that appeared to show that 
unrealistically small optical depths of neutral hydrogen in the ISM were required to allow ionizing photons 
to reach gas at large distances from the O stars. However, the analysis by Seon ( [2009J did not comprise 
a correct photoionization simulation where the ionization state of the gas was calculated self-consistently. 
Instead the penetration of ionizing photons through a neutral ISM was characterized using an exp(— r) factor 



where r is the optical depth due to neutral hydrogen. Seon (2009) correctly stated that the optical depths 



of neutral hydrogen in a smooth ISM are large and will not allow photons to reach large distances from 



the midplane. Instead, Seon (2009) claims that optical depths of neutral hydrogen need to be smaller by 
factors of around 10~^ to allow ionizing photons to penetrate from midplane O stars to high altitudes in 
M51. However, in a photoionization simulation the neutral hydrogen fraction in diffuse ionized gas is in fact 
very low, typically 10~^ and even smaller closer to ionizing sources ( Wood fc Mathis|2004 |. Therefore when 
the medium is ionized Lyman continuum photons can indeed penetrate to large distances. In addition, as 
demonstrated in ours and other simulations clumping of the gas in the ISM produces low density regions 
and this further helps the ionizing photons reach and ionize gas at large distances from the midplane. The 
combination of small neutral fractions in the DIG and low density regions in a 3D ISM will combine to 
produce the low column densities of neutral hydrogen (compared to a smoothly distributed ISM) required 
for ionizing radiation from O stars to reach and ionize the DIG in M51. 



5. Conclusions and future work 

Our main conclusion is that the density structures produced by 3D hydrodynamical simulations of a 
turbulent ISM allow ionizing photons from midplane sources to reach and ionize gas at large altitudes. This 
is in agreement with previous photoionization simulations of analytically produced fractal density structures 
that also indicate that a 3D ISM is required for O stars to ionize high altitude gas. We therefore believe 
that the results of the photoionization simulations of analytic ISM densities and now those of more realistic 
hydrodynamical simulations lend further compelling support for O stars being the dominant ionization source 
of the DIG. 

While a 3D ISM can solve the problem of propagating ionizing photons to large altitudes, our results 
indicate that the dynamical simulations likely require additional physics to reproduce the details of the DIG 
observations within the Galaxy. Gompared to the Galactic DIG, the simulations in this paper show broader 
EM distributions suggesting that magnetic fields may be the missing ingredient to reconcile observations 
with theory. Magnetic fields will reduce the density variations and provide pressure to support gas and give 
a higher density at kiloparsec heights above the midplane. This is borne out in the ISM MHD simulations of 



de Avillez & Breitschwerdt ( 2005 ) and preliminary simulations of our own that incorporate magnetic fields 



(Hill et al. in preparation). 

The hydrodynamical simulations in this paper use the FLASH v2.4 code developed by the Center for 
Thermonuclear Flashes at the University of Chicago (Fryxell et al. 2000). MRJ and M-MML are partly 



supported by NASA/SAO grant TM0-11008X from the Chandra Theory Program. The Wisconsin H- Alpha 
Mapper, ASH, LMH, and RJR are supported by the National Science Foundation through grant AST-06- 
07512. GJM is supported by a Research Fellowship from the University of Sydney 
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Fig. 1. — Mean density (solid line) from the hydrodynamic simulations compared with the Dickey- Lockman 
distribution for neutral hydrogen (dashed line) and a Dickey-Lockman plus extended component representing 
diffuse ionized gas (dot-dash line). 



12 





Fig. 2. — Average neutral fraction as a function of height for ionizing luminosities (left panel) Q49 = 0.1 
(black), 0.5 (red), 1 (green), 2 (blue) and (right panel) Q4g = 2 (black), 3 (red), 4 (green), 5 (blue). 
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Fig. 3. — Edge-on visualizations of the simulation grids. Column densities of neutral hydrogen (top row) 
and ionized hydrogen (middle row), and emission measure {J ri^ds) (bottom row) are shown for ionizing 
luminosities of (left to right) Q49 = 0.5, 1, 2, 4. 
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Fig. 4. — Mean densities of neutral (left panel) and ionized hydrogen (right panel) for ionizing luminosities 
Q49 = 0.5 (black), 1 (red), 2 (green), 4 (blue). Overplotted with dashed lines are the Dickey- Lockman 
distribution for neutral hydrogen (left panel) and the average electron density derived for the Galactic DIG 
(right panel). 
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Fig. 5. — Neutral fraction as a function of height through the center of the grid for simulations with a 50 pc 
spherical blob located at 1.5 kpc above the midplane. The ionizing luminosity is Q4Q = 2s^^ and the curves 
(lowest to highest) are for blob densities of 0.01, 0.05, 0.1, 0.5 cm~^. As the density increases the blob 
transitions from being almost fully ionized, to having an ionized surface skin, to being almost fully neutral. 
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Fig. 6. — Emission measure distributions for the Galactic DIG from the WHAM survey and from our 
(549 = 2 photoionization simulation. Fraction (dex)^^ is the fraction of points in each bin divided by the 
logarithmic bin interval. Left panel: The black histogram is of J n\ds in the simulation, viewed face-on. The 
red histogram is of EM sin |&| from the WHAM survey with classical H II regions and sightlines with |6| < 10° 



removed, from Hill et al. (2008). Lognormal fits to each distribution are also shown. Right panel: The black 



histogram is of J nids in the simulation, viewed edge-on. The red histogram is of EM cos |5| from the WHAM 
survey toward the Persus Arm for velocities in the range —100 km s~ < I^sr < —25 km s^ with classical 
H II regions removed. For the emission measures from the simulations we have removed contributions from 
midplane regions with \z\ < 150 pc. Notice that the model distributions are broader than the WHAM 
observations. 
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Fig. 7. — Histograms of sonic Mach number for subsets of the hydrodynamical model. Black line: entire 
sample. Green line: all gas at \z\ > 125 pc. Red line: all warm (defined as 5000 < T/K < 20000) gas. Blue 
line: warm gas at \z\ > 125 pc. The sound speed in each cell is calculated as v — 7.18 (T/8000 K)^/^ km s^ , 
appropriate for a solar metallicity, neutral gas. 



